First order sidebands in circuit QED using qubit frequency modulation 
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Sideband transitions have been shown to generate controllable interaction between superconduct- 
ing qubits and microwave resonators. Up to now, these transitions have been implemented with 
voltage drives on the qubit or the resonator, with the significant disadvantage that such implemen- 
tations only lead to second-order sideband transitions. Here we propose an approach to achieve 
first-order sideband transitions by relying on controlled oscillations of the qubit frequency using a 
flux-bias line. Not only can first-order transitions be significantly faster, but the same technique can 
be employed to implement other tunable qubit-resonator and qubit-qubit interactions. We discuss 
in detail how such first-order sideband transitions can be used to implement a high fidelity CNOT 
operation between two transmons coupled to the same resonator. 



I. INTRODUCTION 

In circuit QED, superconducting qubits are coupled to 
a common mode of a transmission line resonator [TJ |2] or 
3D cavity [3]. The shared mode can be used to mediate 
interaction between the otherwise non-interacting qubits. 
Several mechanisms to realize two-qubit gates in this sys- 
tem have been proposed [3H3 and implemented [BTI12). 
One approach relies on sideband transitions [4] where, 
for example, the qubit state can be swapped with the 
resonator state. By combining sideband transitions in- 
volving two qubits, an entangling two-qubit gate can be 
implemented [13l [14] . 

In practice, sideband transitions are realized by 
voltage-driving the resonator [H [TS] or the qubit [21 Qj5] . 
Because this type of drive can only cause transitions be- 
tween states of different parity, and since the two states 
involved in any sideband transition have the same par- 
ity, voltage driving can only be used in a second order 
process j2J. In other words, in circuit QED, all sideband 
transitions studied so far have relied on two-photon pro- 
cesses, leading to low transition rates and consequently 
slow two-qubit gates. 

In this paper, we propose to speed up gates based on 
sideband transtions by relying on a different kind of drive. 
Indeed, in many current experiments [8, 17 19 , the qubit 
transition frequency can be adjusted by several hundreds 
of MHz in a few nanoseconds. This is made possible by a 
fast- flux line fabricated in close proximity to the qubit 's 
loop and is described by the Hamiltonian 

Here, f(t) represents the qubit frequency change and is 
a function of the externally applied flux. By modulating 
the function fit) cx cos(wfc^): and choosing the modula- 
tion frequency ojfc appropriately, we show that sideband 



transitions can be generated. Since this frequency control 
(FC) signal is proportional to er z (and not a x as is the 
case for the more typical voltage driving) , sidebands can 
be generated in a first order process and consequently re- 
sult in faster two-qubit gates. Moreover, in the presence 
of two-qubits coupled to the same mode of a resonator, 
this frequency modulation can result in a second-order 
flip-flop interaction between the qubits in a way simi- 
lar to the FLICFORQ [2D] and the cross-resonance [?T] 
gates. 

This paper is divided as follows. In Section [H] we give 
the Hamiltonian for two qubits coupled to the same mode 
of a resonator, including an FC drive on both qubits. By 
moving to the frame that diagonalizes the Hamiltonian, 
we find that the FC drive generates multiple effective 
qubit-resonator and qubit-qubit interactions. These in- 
teractions can be essentially turned on and off by choos- 
ing the appropriate modulation frequency wfc- Here, we 
focus on sideband transitions which are shown to be of 
first order in the modulation amplitude. In Section |III[ 
we calculate the dynamics associated to the red sideband 
process between a qubit and the resonator, and consider 
the error coming from the interaction with a spectator 
qubit. These ideas are then extended in Section |IV| to 
the more realistic case of many-level systems (MLSs) in 
a resonator. We also show that the presence of a second 
excited state can be exploited in a pulse sequence con- 
taining red sidebands to generate a two-qubit gate which 
is locally equivalent to a CNOT. Finally, in Section [V] 
we explain how our proposal can be implemented with 
transmon qubits. We first consider the case of a sin- 
gle red sideband and show how possible caveats can be 
avoided. At last, wc present simulation results that pre- 
dict that FC modulations can be used to generate high- 
fidelity maximally entangling gates with two transmons. 

II. EFFECTIVE TUNABLE INTERACTIONS 
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In this section, we first present the lab frame Hamil- 
tonian for two qubits coupled to a single mode of a res- 
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onator. Assuming the qubits to be well detuned from 
the resonator and ignoring the presence of any driving 
terms for the moment, we eliminate the resonant qubit- 
resonator interaction using a dispersive transformation. 
This results in, in particular, a residual qubit-qubit vir- 
tual interaction which we diagonalize. Reintroducing the 
FC drive Hamiltonians in that diagonal frame, we ob- 
tain the expressions for the tunable interactions, among 
which the red and blue sideband transitions are present. 
To simplify the discussion, we will ignore here the pres- 
ence of higher excited states of the qubits. These levels 
are important, for example, in the transmon qubit and 
their effect will be discussed in Section Hvl 



is a qubit-qubit interaction mediated by the exchange of 
virtual photons. 



B. Diagonal Hamiltonian 

Because of this last term, ifo is not diagonal. It how- 
ever has a simple structure which we now exploit. Indeed, 
since Hd is diagonal with respect to the cavity degree of 
freedom, it is useful to focus on subspaces labelled by the 
fixed photon number n. In these subspaces, the Hamil- 
tonian effectively reduces to a 4 x 4 matrix 



A. Lab frame Hamiltonian and generalized 
dispersive Hamiltonian 

The electric-dipole coupling of two qubits to a single 
mode of a resonator is described by the Rabi Hamilto- 
nian, with h = 1, [T] 



with 
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H = uj r a) a + 



k=l,2 



,(*) 
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Aq = A q 



-2nA s 
2n£ s 



(8) 
(9) 



Here, uj r is the resonator frequency, u>a the transition 
frequency of qubit k, <?( fc ) th e qubit-resonator coupling 



strength for qubit k and X 



,t 



a. The free qubit 



eigenstates will be labelled {\g}, |e)}. 

To approximately diagonalize H we use the generalized 
dispersive transformation based on the unitary [351 123] 



J7d = exp 



A (fc W + A^W +^a k a 2 -B.c. 



fe=l,2 



,(3) 



where A<*> = cjW/AW, A< fe ) = .g( fe )/£( fc ) with = 

uja — w r and = &a + w r- The squeezing-like 

terms are proportional to £}- k ^ = (x^ +M < '' c ' ) )/4w r , where 
^(fc) _ g(k)^(k) an( j ^(fc) _ g(fe)A(/c) are respectively the 

dispersive and Bloch-Siegert shifts. Using the Campbell- 
Bakcr-Hausdorff relation 



-x i 



#e A =H+[H, X] 



[[H,X],X] 



(4) 



and where we have defined Aq = Cba — uj a , Ag = 
- ^(2) , Sq = qO) + ftCO ? and Ss = 5(1) + 5(2) . B e- 
cause of its block structure, this Hamiltonian can easily 
be diagonalized. As in the cross-resonance gate [61 |2"T]. 
we will define two effective (or logical) qubits out of the 
eigenstates of this 4-dimensional block, the dependence 
on n simply leading to ac-Stark shifts of these effective 
qubits. As discussed in details in Appendix [AJ if the 
two qubits are far-detuned from each other, the resulting 
diagonal Hamiltonian is, to second order in J/Aq and 
J/^q, given by 

tfdiag = u r (Ja + (u>P + 2a t aS (1) + S+) ^ 

+ {^+2a^aS^+Sj) T l (10) 

where 



we obtain, to second order in the couplings g^ k \ 



UlHUu 



~(fc) 



w r aW ^ r <T k z +S^a k z a^a + Ja 1 x al (5) 



fe=l,2 



where u)i 



,(fe) 



t^a ' +X^ fe ) +^i' fe ' are the Lamb-shifted qubit 
frequencies and S^ k ' — x + M is the ac-Stark shift 
per photon. The last term of Hy>, proportional to 



J = 



? (i) 9 (2) 



fc=l,2 



1 



1 



A( fc ) E( fc ) 



(6) 



5' 



(ii) 



Here, we have also defined the Pauli operators T k act- 
ing on the effective qubits. In this diagonal frame, the 
interacting physical qubits are replaced by two logical 
non-interacting (and non-local) qubits, whose eigenstates 
will be denoted {|0), |1)}. The effect of the resonator- 
mediated J interaction between the physical qubits is to 
shift the logical qubits by a quantity S'j which is as- 
sumed here to be small (see Appendix |A| for the general 
result). Since it is of fourth order in 5™, this shift can 
be neglected, as will be done in the rest of the paper. 
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C. FC drive in the diagonal frame 

We now introduce an FC drive on each qubit. In the 
lab frame, this takes the form 



H FC (t) 



f {1) (t) t , 



(2)/ 



Later on, we will take 

y(fe)( t ) =e (*)( t ) C0B 



(12) 



(13) 



where ' (t) is an envelope that varies slowly compared 

(k) 

to the modulation frequency Wpg. Going to the disper- 
sive frame and then to the diagonal frame as above, this 
Hamiltonian becomes 



+ H PO (t)+H QQ (t). 



(14) 



The first two terms represent modulation in the effective 
qubit splittings 



where we have defined 



2 

z£ 

2 



s(fe) 



At = 



(A( fe )) 2 + (A< fc )) ; 
J 



(a f a + 1/2) , 



A Q + 2ataA s 



(15) 
(16) 

(17) 
(18) 



As will be seen in Section VB[ these terms can lead to lo- 
cal phase gates on the logical qubit. Just like in the single 
qubit case 23J, in addition to these frequency modula- 
tions, we also get sideband transitions and parametric 
oscillations. However, these processes are modified by 
the qubit-qubit interactions. First, we find 



A«aV +A< 1 Vi 



+ AJT* (\WJt 2 _-AWJtI) +H.c. 



A< 2 Vr 2 -f- A^aM 



(19) 



(20) 



A. 7 (A( 2 >aV +A( 2 ) a Vi) 



H.c. 



corresponding to red and blue sidebands. In practice, 
these terms will average out unless the modulation fre- 
quency wfc is chosen appropriately. For example, tak- 
ing — A( fc ) or ~ £( fc ), respectively lead to 
red or blue sideband transitions between the resonator 
and qubit k. These transitions respectively correspond 
to |1; n) <-> |0; n + 1) and |l;n+l) <-> |0; n), with the first 
number labeling the state of qubit k and the second be- 
ing a Fock state of the resonator. Additionally, because 



of the qubit-qubit interaction J, we also get leakage-like 
terms which are proportional to Aj: driving qubit k at 
the frequency corresponding to a sideband for qubit k' 
will drive the sideband in qubit k' . In this process, qubit 
k will pick up a phase. While these terms could be useful 
for logical gates, they are in practice very small because 
of the Aj factor which is of second order in g( k \ Fur- 
thermore, the assumption of a large qubit-qubit detun- 
ing makes these leakage terms far off-resonant with the 
other sideband terms, making these spurious transitions 
negligible in practice. We also point out that in most 
cases, A' fe ) ^> A^ k \ implying that the red sidebands rate 
is much larger than for the blue sideband. 

As in the single qubit case, we also unsurprisingly 
find a term which corresponds to parametric amplifica- 
tion [23] 

ffpo(i) -~ X (k) A {k) f w (t) r k (a 2 +a t2 ). (21) 

fe=l,2 

This term is only relevant for cjpc — 2uj r . However, in 
typical situations, it is not large enough to be useful as 
a practical source of squeezed microwave light Q 

Finally, the virtual qubit-qubit interaction is responsi- 
ble for the appearance of a new term in the FC Hamil- 
tonian in the diagonal frame 



H, 



3Q( ^{_/ (1) ft)+/ (2) (0 ( A ( 1 ) A ( 2 )_ A ( 1 ) A ( 2) ) 

-A,,[/W(i)-/ (2) W]}(rlr 2 +r|r!) (22) 
/(1)( ^ /(2)W (A« A< 2 > - A(i>A( 2 >) (rW + r\rl 



Because of this term, when wpc — Aq, a flip-flop inter- 
action between the two logical qubits is turned on (first 
two lines of Hqq). This can be used to generate VSWAP 
gates to entangle the qubits. Special care must however 
be taken when it comes to the relative phase between the 
modulation functions. Indeed, taking f^\t) and f^ 2 \t) 
in phase results in a cancelation of the term proportional 
to f^'(t) — f( 2 \t), while taking them to be out of phase 
cancels the term in f^'(t) + f^ 2 '(t). Therefore, depend- 
ing on the qubit and resonator frequencies, one should 
choose the phase that allows to select the fastest popula- 
tion transfer. We also point out that if uJa < uj r < 0Ja , 
the signs of X^ A^ 2 ' and A' 1 ) A' 2 ) are opposite, allowing a 
larger Rabi frequency if A( fc )'s are comparable to A^'s. 
We can also drive a transition that creates or destroy a 
quantum in each qubit when wfc — (last term line of 



The standard deviation in the X quadrature of a degenerate 



4e , where 

4 



parametric oscillator when ajpc — 2ui r is AX 
r = |arctanh(4e(' c ' A( fc 'A( fc '/^r) [H]. For the generalized disper- 
sive transform to be valid, we require A( fc ), A( fc ) < 1. Therefore, 
r > 1 implies e' fc ' S> w r , which is impossible if ui a ~ w r as is 
usually the case in circuit QED. 
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Hqq). Because of its g 4 dependance, the effect of Hqq 
is small for the parameters chosen in this paper and will 
be ignored in the following. 



III. SIDEBAND CONTROL WITH FC DRIVES 

Since they are the terms of largest amplitudes in Hp^ s 
that can lead to two-qubit gates, we focus here on the 
red sideband terms present in the first lines of H^ B and 
i?| B . As already pointed out, with respect to previous 
proposals [4] and realizations [Mj [16] in circuit QED, 
we emphasize that the sideband rate under an FC drive 
is of first order in <?, rather than of second order. In 
this section, we calculate the evolution of the state of a 
target qubit k and of the resonator under a red sideband 
transition generated by an FC drive on k, the other qubit 
kl merely acting as a spectator. However, because of 
the qubit-induced frequency pull of the resonator, the 
resonance frequency for a sideband on the target qubit 
k depends (weakly) on the state of the spectator qubit 

k' 



kl . This leads to a small error that increases with A 
which is evaluated below. Finally, because our goal is to 
see how good a qubit-qubit entangling gate based on the 
red sideband could be in principle, we will ignore here 
the effect of damping and dephasing. These will simply 
lower the fidelities in a straightforward manner. 

As discussed in Section [TXJ at the appropriate FC fre- 
quency, red sideband transitions occur between eigen- 
states in the diagonal dispersive frame. These eigen- 
states can be written in the basis {|tiT2;ti)}, where 
r k e {0, 1} denotes the eigenstates of logical qubit 
k and n represents Fock states of the dispersively- 
shifted resonator. Focussing on red sideband tran- 
sitions, which preserve the total number of quanta, 
the Hilbert space can be restricted to the subspace 
{|00;n + l),|10;n),|01;n + l),|ll;ra)}. To simplify the 
discussion, we keep only the relevant terms in the total 
transformed Hamiitonian in the presence of an FC drive 
at the red sideband frequency. This leads to 



wi k) + 2at ffl S< fc > 



-(fc) 



A«/W(t)(atrl 1) +ar| 1) ). (23) 



This is a good approximation when A^e' 1 - 1 <C \u)^ — 
oj r \, \uj^ — w r \, and \&a — <Da |, i.e. the transitions that 
can be stimulated by the FC drive are well separated. 
Furthermore, we have neglected Hl(t) and H%(t) which 
correspond to oscillations in the qubit frequencies on a 
timescale that is much smaller than the sideband transi- 
tions. These terms can be safely dropped if the average 
of f^(t) during a whole period is zero, which is the case 
here with a drive described by Eq. (13). 



spect to 
U TO t{t) = exp 



■(*) 



(24) 



(25) 



The resulting Hamiitonian is 

H' cd (t) = -A«/W (t) [ar^ exp (zA«i) + H.c. 

where A^ 1 ] = ^ - uj r + 2S^a^a - S^tP. Working 
in the restricted Hilbert space and ignoring dissipation, 
Schrodinger's equation for the above Hamiitonian and 
with the drive defined by Eq. (13), leads to two uncou- 



pled sets of differential equations for the probably am- 
plitudes c^i x {t) and corresponding to the second qubit 
being excited (+) or not (— ) 



Cq (t) = icf(t)e n 
cf(t) = icg(t)e n 



e i(«po-A±)i , e -i(w FC +A±)i 
e -i(o) FC -A±)i + e i(w FC +A±)t 



where wpc = anc ^ 



(26) 
(27) 

(28) 
(29) 



For e„ -C u) + A* , the fast oscillating terms can be ne- 
glected and an analytic solution can be found. The evo- 
lution operator corresponding to this solution is, for each 
subspace ± 



v ± (t) 




cos ■ 



■ <5± • rt 

l T 81n I 



e l5 ±*/ 2 sin f 



-2i 

e iS±t/2 



iS±t/2 



sm ■ 



cos 4r + i— sm 4f 



(30) 



where r = y 6± + 4e 2 and 5± — wpc — A^ is the detun- 
ing between the FC drive and the Stark-shifted resonance 
frequency of the red sideband. We already notice that 
without coupling between the second qubit and the res- 
onator, we eliminate the Stark shift and thus can have 
<5 + =5- =0. This leads to full population inversion 
between |0;n + 1) and |l;n) at the Rabi frequency 2e n , 
which is of first order in g\ k ) / /\S- k > . 

Given the evolution operator V ± (t) , we now calculate 
the average fidelity to a red sideband 7r-pulse on the first 
qubit, irrespective of the state of the second qubit. The 
ideal process is described by 



It is useful to move to the interaction picture with re- 



U 



-»(|Q;n + l)<l;n| + |1; n)(0; n + 1|) ® I Q2 . (31) 
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To compute the average fidelity, it is useful to define the 
Choi matrices corresponding respectively to the ideal Ctj 
and the actual Cy(t) processes [2"5"l |2"6] 

Cu = \ £ \W\c®U\t)(j\tf (32) 

i ,je{a,b} 

Cv{t) = \ J2 \i){3\c®V(t)\i){j\VHt). (33) 

i,je{a,b} 

The sums above are evaluated over states \a) — |0; 1) 
and | b) = |1;0) that are involved in the red sideband 
quantum process. The subsystem labelled c in the above 
equations represents a copy of the relevant Hilbert space 
{|0; 1), |1; 0)} and is required in the definition of the Choi 
matrix. The gate fidelity Fuv(*) = tr [C\jCy(t)] of that 
process is thus 



1 



Fuv(«) = 7 |tr [C/V(t)]| 



Using Eq. (32) and (33), we obtain 



rt 



Fuv(t) = 2-| sin 2 - (1 + cosM) , 
r l 2 

from which we obtain the average fidelity [37] 



Fuv — 



2F UV + 1 



(34) 



(35) 



(36) 



As expected, the fidelity is unity with e n t = tt/2 and if 
the Stark shift S^ 2 ' due to the second qubit vanishes. In 
practice, reduces the average fidelity. 

As an example, we now assume the modulation fre- 
quency to be at the red sideband transition frequency, 
given that the second qubit is in its ground state. In 
other words, we have 



. if the second qubit is in state |0) , 
± = i 2S^ if the second qubit is in state |1) ' ^ ' 



In this situation and for a given temporal shape of the 
FC modulation, the condition for population inversion 
becomes e n t = 7r/2, with 



1 



i(*)dt. 



(38) 



p Jo 



Taking e n t p = n/2 and replacing e„ — > e n in Eq. (35) 



yields a simple expression for the gate fidelity in the limit 

S+ < e„ 



2ei 



4ef 



(1 + cos S+tp) . 



(39) 



In words, the infidelity 1 — Fuv is minimized when the 
Rabi frequency that corresponds to the FC drive is large 
compared to the Stark shift associated to the spectator 
qubit. The average fidelity corresponding to the gate fi- 
delity Eq. ( 39 1 is illustrated in Fig. [I] as as a function 




3 4 5 6 
S*( 2) (MHz) 



FIG. 1. (Color online) Average error with respect to the per- 
fect red sideband process |1; 0} o j0;l). A gaussian FC pulse 
is sent on the first qubit at the red sideband frequency as- 
suming the second qubit is in its ground state. Full red line: 



average error of the red sideband as given by Eq. ( 39 1 when 
the second qubit is excited. Blue dashed line: population 
transfer error 1 — Pt, with P t given by Eq. (41 1. Black dots: 
numerical results for the average error. We find the evolu- 
tion operator after time t p for each eigenstate of the second 
qubit. The fidelity is extracted by injecting these unitaries in 
Eq. (34). The qubits are taken to be transmons, which are 



modelled as 4-level Duffing oscillators (see Section VAl with 



E { j } = 25 GHz, Ej 
MHz, yielding 
J 1 ) 



(2) 



35 GHz, E ( c ] = 250 MHz, E% 



(2) 



5.670 GHz and ui, 



(2) 



300 

7.379 GHz, and 

<7oi' = 100 MHz. The resonator is modeled as a 5-level trun- 
cated harmonic oscillator with frequency ui r — 7.8 GHz. As 
explained in Section [V] the splitting between the first two lev- 
els of a transmon is modulated using a time-varying external 
flux (j>. Here, we use gaussian pulses in that flux, as described 
by Eq. (71 1 with r = 2a, a — 6.6873ns, and flux drive am- 
plitude A(p — 0.075 00- The length of the pulse is chosen to 
maximize the population transfer. 



of 5( 2 ) (red line) assuming the second qubit to be in 
its excited state. We also represent as black dots a nu- 
merical estimate of the error coming from the spectator 
qubit 's Stark shift. The latter is calculated with Eqs. ( 34 ) 



and (36). Numerically solving the system's Schrodinger 



equation allows us to extract the unitary evolution op- 
erator that corresponds to the applied gate. Taking U 
to be that evolution operator for the spectator qubit in 
state |0) and V the operator in state |1), we obtain the 
error caused by the Stark shift shown in Fig. [T] The nu- 
merical results closely follow the analytical predictions, 
even for relatively large dispersive shifts S'' 2 -'. 

Reducing the dispersive shift of the spectator qubit 
S<- 2 > rapidly increases the average fidelity. There are sev- 
eral ways to do this. The straightforward solution is to 
bring it to a very large detuning with the resonator, thus 
reducing the Stark shift. One could also use a tunable 
coupling [23 [33] . Finally, if neither of these options is 
convenient, this kind of issue can typically be addressed 
by pulse shaping techniques [30 . 

It should also be noted that if one is interested only 
in population transfers and not in phase information, 
the above analysis exaggerates the error because it cares 



G 



about phase. In this context, a more appropriate figure 
of merit is the population of the target state 



P t = \(l;n\V(t)\0;n + l)\' 



4r 2 



(40) 



Taking, as in the previous case, t = t } 
expanding to second order in 5+/e n , this reduces to 



7r/2e„ and 



4r- 



5\ 



Ac 



■2 ' 



(41) 



This expression is the blue dotted line in Fig. [T] Clearly, 
not taking phases into account yields a much smaller er- 
ror rate. As will be discussed in Section [V C| this feature 
remarkably helps the implementation of two-qubit gates, 
Indeed, in that case, the specific phase of each sideband 
gate is irrelevant, provided that it is well-controlled and 
stable over many repetitions of the experiment. 



IV. MANY-LEVEL SYSTEMS AND LOGICAL 
GATES BASED ON SIDEBAND CONTROL 

While the two-level model presented in the previous 
sections is useful because of its simplicity, it is not an 
appropriate description of most superconducting qubits, 
such as the transmon [31 or the CSFQ [32]. Indeed, 
these devices are best described as many-level systems. 
In this section, we generalize the previous results to that 
more realistic situation and discuss how these additional 
levels can by exploited to our advantage. As above, we 
start be obtaining an effective Hamiltonian using a dis- 
persive transformation and then diagonalize the remain- 
ing qubit-qubit interaction. We only outline the main 
steps of this calculation, the details of which can be found 
in Appendix [B] before presenting a sequence of pulses 
corresponding to a qubit-qubit entangling gate that takes 
advantage of the many-level structure of the qubits. 



A. FC drive Hamiltonian for many-level systems 

The Hamiltonian of two many-level systems (MLS) 
coupled to the same harmonic mode takes the form [3T] 



H = w r cta + Y, E w I fe)nW 



I 1,1 



(42) 



fc=l,2 ieHt 



+ E E ^(n^U + n^^+H-c. 



with n 



fc=l,2 ieV.M-1 



(fc) 



,0) 



\i)(j\ {k) acting on the fc-th MLS, u\ 

(k) 

the bare MLS frequencies, g\ ' the coupling rate be- 
tween states i and i + 1 of the fc-th MLS and the 
resonator, and M the number of levels in each MLS. 
T~Lm = {Iff), |e), |/), | ft), ...} is the Hilbert space of a bare 
M- level system. For simplicity, we have assumed that 



only nearest-neighbor states are coupled by the electric- 
dipole interaction. As will be discussed further below, 
this approximation is well satisfied for transmons 3 1 J . 

As shown in Appendix [B] the dispersive transforma- 
tion described by Eq. ([3| can be generalized to the multi- 
level case. This transformation f/|f LS , given by Eq. (Bl I, 
is used to approximately diagonalize the above Hamilto- 
nian, leading to Stark shifts. As in the simple case, this 
also leads to a qubit-qubit coupling term of amplitude 
J. An additional transformation [/MLS is introduced to 
diagonalize this effective coupling. Successively applying 
U^ LS and Uf LS , we find to second order in gf 1 



diag _ 
D — 



oj r a) a + 



M-l 

E E 

k=l,2 i=0 



o( k )-n( k ) t 



where 



-(fc) (k) . r (fc) 
7-0) (k) (fc) 

c (fe) (k) . (k) (fc) (k) 



(43) 

(44) 
(45) 
(46) 



Here, we have defined x 

,(*)* A (k) 



(k) 



(k)*\(k) j (k) 

91 K and fj,\ 



g\ ' A\ , respectively the dispersive and Bloch-Siegert 
shifts with Af> = ft « /A| fe », Af> = fff/sf, = 



,(fc) , ,(fc) 



lu^—lu- ' —uj r , and s| fe = —wl"' +uj r . The eigen- 
states of the above Hamiltonian are, in each photon sub- 
space, logical (dressed) transmon states |7ri7T2;n), where 

ir k e {0,1,2,...}. 

We now introduce the FC drive on each MLS. The 
most general Hamiltonian that describes this in the bare 
frame is 



,(fe) 



,(*) 



-= E E f? ] w%, 

fc=l,2 ieWj; 



(47) 



where fi k \t) are real functions of time. Successively 
applying U^ hS and Uj 11 ^ 5 on that Hamiltonian, we ob- 
tain an expression for Hpc in the frame that diagonal- 
izes Eq. (42). To simplify the discussion, we drop here 



ffc) 

terms that are of second order in g\ because they lead 
to tractable but unsightly expressions. In this way, we 
find that the first-order corrections are sideband terms 



ti 



diag 
FC 



fc=l,2 



M-2 

E 

=0 



(fc) 

SB,i' 



(48) 



where 



4% = ^ k \t) [Xf^l^ + Af )Hg +1 a + H.c. 

(49) 

and 5f?\t) = f}»(t)-fV ,) (t). 

Thus, instead of having only one red and one blue 
sideband transitions for each photon number as in the 



7 



-|2> 




vr = o) K = i) k = 2) 



FIG. 2. (Color online) Sideband transitions for a three-level 
system coupled to a resonator. Applying an FC drive at fre- 
quency Ai,j+i generates a red sideband transitions between 
states \i + l;n) and \i;n+ 1), where the numbers represent 
respectively the MLS and resonator states. Similarly, driving 
at frequency E iji+ i leads to a blue sideband transition, i.e. 
\i; n) o |i + l;n+l). Transitions between states higher in 
the Fock space are not shown for reasons of readability. This 
picture is easily generalized to an arbitrary number of levels. 



two-level case, here M — 1 red and blue sidebands are 
possible for each n. These are illustrated in Fig. [2j For 
large enough anharmonicity of the MLSs, these sidebands 
are at distinct frequencies. As will be discussed in the 
next section, this selectivity is a useful resource to design 
quantum logical gates. 



B. Pulse sequence generating a maximally 
entangling gate 

As already mentioned, sideband transitions can be 
used to generate entanglement between two qubits cou- 
pled to the same resonator. Here, we show how the sec- 
ond excited dressed level of a MLS, labelled |2), can be 
exploited to design a maximally entangling gate. 



Let us define a 



n 



n 



(k) 



. - hJ , ^ , corresponding to 

a 7r-pulse between the levels i and j of the k-th logical 



MLS, and R 



(fc) 



l {k) P t 



n (fc) a 

11 i+l,i u ' 



corresponding to 



a full red sideband transition between the fc-th MLS and 
the resonator. Using these definitions, we introduce the 
pulse sequence 



77 _ d(i)d(2)_(2) r( 2 )rW 

Cent — -n-oi -""12 L, £c,12 rl 12 -"01 • 



(50) 



Up to phases, this is equivalent to the CNOT. Indeed, 
looking at the effects of this sequence on each state of 
the two-qubit computational basis {|00), |10), |01), |11)} 
given that the resonator is initially in its ground state 
and ignoring phases for the moment, we have 

|00;0) -> |00;0) -> |00;0) -> |00;0) -> |00;0) -> |00;0) 
|10; 0) -> 1 00; 1) -> |00; 1) -> |00; 1) -> |00; 1) -> 1 10; 0) 
|Q1;0) |01;0) |01;0) ->■ |02;0) |01; 1) -> |11;0) 
|11; 0) -> |01;1) -> |02;0) -J- |01;0) -> |01;0) -> |01;0). 



Therefore, the first qubit is flipped if the second one is 
in state 1 1) . This is represented more generally in matrix 
form by 



L'fint, — 







i<t>i 



(I 




\0 



\ 



e^ 2 

g»03 o J 



(51) 



where fa, fa, and fa are phases that depend on the shape 
of the FC pulses. Indeed, as will be seen in Section [VB| in 
realistic situations, the FC drive is not strictly a cosine 
drive on all the qubit transitions frequencies, but also 
has a component that can be interpreted as a constant 
shift on these frequencies, leading to additional phase 
accumulations. 

The above unitary is equivalent, up to one-qubit gates, 
to the CNOT. Indeed, we can write 



CNOT = u£ ) U eat U™u£\ 



(52) 



(k) 

where Ug is a rotation of qubit k along the z-axis by an 
angle 6 and 



9i = (fa -fa- fa)/2, 
% = (fa -fa- fa)/2, 
h - (fa -fa- fa)/2. 



(53) 
(54) 
(55) 



It should be noted that the above sequence should be 
faster than the one presented in Ref. [I], since it con- 
tains four sideband pulses instead of five. Together with 
the fact that sidebands generated with FC drives have 
a first-order nature, this allows much faster results than 
expected in the initial proposal. 



V. PHYSICAL IMPLEMENTATION WITH 
TRANSMONS 

In the previous sections, we first introduced FC drives 
on two-level systems, and then of generic MLSs. Here, 
we focus on a specific kind of MLS: the transmon 
qubit [3TJ [33]. Using a simple model, we predict that 
the red sideband transition can be generated with a very 
high Rabi frequency using an FC drive on that qubit. 
We also show simulation results for the pulse sequence 
presented in Section [IV B| 



A. The transmon as a Duffing oscillator 

A simple Hamiltonian for the transmon qubit can be 
derived from the Cooper-pair box Hamiltonian by ex- 
panding it to fourth order in Ec/Ej. In the limit where 
Ec/Ej <C 1, the system is well approximated by a Duff- 
ing oscillator [3T] 



H t , 



^/SEcEjb^b 



Eg 
12 



(56) 



where Ej and Ec are respectively the Josephson and 
capacitive energy of the transmon, while b and W are 
the ladder operators of this weakly anharmonic oscillator. 
The associated cigcnfrcqucncics are then 



Ec\ . E c . 2 



(57) 



with uip — \/8EcEj the plasma frequency which depends 
on Ej. Additionally, the transmon- resonator coupling 
rates are, again in the regime of large Ej / Ec ratios, given 
by 



9j,j+l - 9geVT+l, 



where 



-2i/3eK°m 



Ej 
V2 V8^< 



(58) 
(59) 



(60) 



In addition, gj+i,, — g*a + i- The transmon-resonator cou- 
plings also depend weakly on Ej. 



B. FC drives on the transmon-resonator system 

If the transmon's Josephson junction is replaced by a 
SQUID, Ej can be tuned with an external flux. Assum- 
ing the flux modulation to have the form of a cosine wave, 
we have 



Ej(t) = Ejx cos [fa + A0cos(o;fc*)] , 



(61) 



where fa and A(f> are the mean value and amplitude of 
the external flux drive. All the fluxes are expressed in 
units of $o/7r, where 4> is the flux quantum. This re- 
sults in modulations of the plasma frequency and thus in 
the transmon energies, something which can be used for 
FC driving. Additionally, this also leads to modulation 
of the transmon-resonator couplings. These modulations 
qualitatively have the same effect as the transmon fre- 
quency modulation, leading to red and blue sidebands. 
However, because of the weak dependence of the cou- 
pling on Ej, this is a small effect which can safely be 
dropped. 

Having dropped the modulations of the coupling, we 
now focus on better understanding the frequency mod- 
ulations. Expanding Eq. (61) around fa to fourth order 
in A<f> and for fa such that tan^i < 1, i.e. fa < 7r/4, we 
obtain the following formula for the transmon transition 
frequencies 



We have also defined u>' p — y/8EcEjs cos fa, the plasma 
frequency associated to the operating point fa. This fre- 
quency is illustrated by the black dots for two operating 
points on Fig. [3^). In addition, there is a frequency shift 
G, standing for geometric, that depends on the shape 
of the transmon energy bands. As is also illustrated on 
Fig. [3^,) , this frequency shift comes from the fact that the 
relation between 0Jjj+i and <fi is nonlinear, such that the 
mean value of the transmon frequency during flux mod- 
ulation is not its value for the mean flux fa. To fourth 
order in A<f>, it is 



G 



1 



tan 2 fa \ , A<j>- 



(4 + 20 tan 2 fa + lb tan 4 fa) oj 



, A^ 4 
1024' 



(64) 



Finally, we get in Eq. ( 62 1 a modulation at each multiple 



mwFc of the flux drive frequency. The amplitude of these 
modulations diminishes with increasing m as 



A</>- 



1 



tan 



A0 3 
16 



tan < 



£3u 



£4u 



tan 2 fa \ , A(j>- 



2 

20 tan 2 



tan^ 



uj' p tan 



15 tan 4 
9% 32 



fa) w 



p 768 : 



(65) 

(66) 
(67) 

(68) 

As illustrated in Fig. [3^,), in the special case where fa = 0, 
i.e. around the flux sweet spot, only the even harmonics 
2cod and Auid are non-zero. Furthermore, in realistic sit- 
uations, A(f> <C 0.1, meaning that e mu falls off extremely 
quickly with m. This allows to focus only on the first 
non- vanishing harmonic coming from the FC drive in the 
following discussion. 

Let m be that dominant harmonic. Then, shaking the 
flux at a frequency such that 



(4 + 20 tan 2 fa + 15 tan 4 fa) u. 



,Afa 
3072' 



mwpc 



A" 



G, 



(69) 



it is possible to induce red sideband transitions between 
states \j + 1; n) O \j;n + 1) at the Rabi frequency 



9j,j+i 



(70) 



4 

£< 

m— 1 



, cos muiFct, (62) 



where Wjj+i = — ojj is the transition frequency 

without the external modulation 



uL-Ec(j + l). 



(63) 



where A" j+1 



to 



3+1 



- ujj —uj r + n(S 



3+1 



Sj is the 



Lamb and Stark-shifted detuning between the resonator 
frequency and the j-th transmon splitting. 

In Fig. J3j}), the Rabi frequencies predicted by the 
above formula are compared to numerical simulations us- 
ing the full Hamiltonian Eq. (42 ), along with a cosine flux 
drive. The geometric shifts described by Eq. (64) are 
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FIG. 3. (Color online) FC driving of a transmon with an 
external flux. The transmon is modelled using the first four 
levels of the Hamiltonian given by Eq. ( 56 1 , using parameters 
Ej/2n = 25 GHz and E c /2tt = 250 MHz. We also have 
g ge /2n = 100 MHz and u r /2n = 7.8 GHz, which translates 
to A ge /27r ~ 2.1 GHz. a) Frequency of the transition to the 
first excited state obtained by numerical diagonalization of 
Eq. (561. As obtained from Eqs. (651 to (68 1, the major com- 



ponent in the spectrum of u ge (t) when shaking the flux away 
from the flux sweet spot at frequency o>fc also has frequency 
ojpc- However, when shaking around the sweet spot, the dom- 
inant harmonic has frequency 2ufc- Furthermore, the mean 
value of w 9e is shifted by G. b) Rabi frequency of the red 
sideband transition |1;0) «-» |0; 1). The system is initially in 
|1;0) and evolves under the Hamiltonian given by Eq. (42 1 
and a flux drive described by Eq. (61 1. Full red line: analyti- 
cal results from Eq. ( |70[ ) with m = 1 and <f>i — 0.25. Dotted 
blue line: m — 2 and <j>i = 0. Black dots and triangles: ex- 
act numerical results, c) Geometric shift for cj>i = 0.25 (full 
red line) and (dotted blue line), d) Increase in the Rabi 
frequency for higher coupling strengths with 4>i = 0.25 and 
A(f> = 0.075. e) Behavior of the resonance frequency for the 
flux drive. As long as the dispersive approximation holds 
(g < g CI it/2ir — 1061 MHz), it remains well approximated by 



Eq. (691, as shown by the full red line. The same conclusion 



also plotted in Fig. 3j:), along with numerical results. In 
both cases, the scaling with respect to A(f> follows very 
well the numerical predictions, allowing us to conclude 
that our simple analytical model accurately synthesizes 
the physics occurring in the full Hamiltonian. It should 
be noted that, contrary to intuition, the geometric shift 
is roughly the same at and away from the sweet spot. 
This is simply due to the fact that the band curvature 
does not change much between the two operation points. 
However, as expected from Eqs. (65) to (68), the Rabi 



frequencies are much larger for the same drive ampli- 
tude when the transmon is on average away from its flux 
sweet spot. In that regime, large Rabi frequencies ~ 30- 
40 MHz can be attained, which is well above dephasing 
rates in actual circuit QED systems, especially in the 
3D cavity [3J. However, the available power that can be 
sent to the flux line might be limited in the lab, putting 
an upper bound on achievable rates. Furthermore, at 
those rates, fast rotating terms such as the ones dropped 
between Eq. (29) and (30) start to play a role, adding 



spurious oscillations in the Rabi oscillations that reduce 
the fidelity. These additional oscillations have been seen 
to be especially large for big relevant e m4) /A™ - +1 ratios, 
i.e. when the qubit spends a significant amount of time 
close to resonance with the resonator and the dispersive 
approximation breaks down. 

Moreover, the analytical model developed here relies 
on the dispersive approximation, and thus breaks down 



unless A 



(k) 



is small. This breakdown is illustrated in 



Fig. [3ji) and [3^), which respectively illustrate the behav- 
ior of the Rabi frequency and the geometric shift with 
increasing g( k '. The regime of validity of the dispersive 
approximation can be captured semi-quantitatively us- 
ing the critical photon number n CT n — A 2 /4g 2 Q]. In- 
verting that relation, we get a critical coupling strength 
g cv it — A/2y / n,. For a red sideband involving only one 
excitation, this reduces to g cr it = A/2. The region for 
which g > g CT - lt corresponds to the shaded green areas in 
the figure. It correctly corresponds to the point where 
our analytical model starts to deviate from the the full 
numerical results. 



C. CNOT gate using FC driving on transmons 

In this Section, we present simulation results for the 
pulse sequence described in Eq. ( 50 ) using the red side- 
band described above for the transmon. The envelope 
function of the FC flux modulation is chosen to be a 
truncated Gaussian 



A<t>(t) = 



^ e -(t-M) 2 /2^ 2 

otherwise. 



if ji — t <t < fi ■ 



(71) 



holds for the Rabi frequency. 



In order to reduce leakage to unwanted levels due to the 
weak anharmonicity of the transmon, we use first-order 

(2) (2) 

DRAG corrections on R\ 2 ' and a x [ 2 pulses [33]. More 
precisely, for each corrected Gaussian pulse, we add a 
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FIG. 4. Amplitude of the gaussian pulse over time. A<j>' is 
such that the areas A+ and 2A~ are equal. Then, driving the 
sideband at its resonance frequency for the geometric shift 
that corresponds to the flux drive amplitude Acj>' allows pop- 
ulation inversion. 



Pulse 


A[<j) /-K,2-K x MHz] 


mH 


u[ns] 


2r[ns] 


R (i) 

R (2) 
R (l) 


0.07308 
0.02520 
51.1412 
0.02520 
0.07308 


16 
46.5 
65.96 
85.42 
115.92 


7 
6.25 
1.48 
6.25 

7 


28 
25 
5.92 
25 
28 



TABLE I. Pulse sequence used in the simulations of a CNOT 
using transmons. The amplitude of the pulses are in units 
of 4>o/ir for the flux pulses and in MHz for the direct drive 
on the second transmon. The total duration of the se- 
129.92ns. Transmon parameters are 



quence is 

25GHz, E ( c ) /2tt = 250MHz, 
= 300MHz, cj) (2) 



AW 



E^/2n 



0.25, Ef ) /2n = 61GHz, 



0.25, giV/2n = g$/2n 



100MHz. Four levels are considered for each transmon. The 
resonator is truncated to 5 levels and has oj r /2n — 7.8GHz. 



drive component that is phase-shifted by it/ 2 and that is 
proportional to the derivative of Eq. ( 71 ) for \t — fi\ < r. 

In practice, using a device such as a transmon, we do 
not only have leakage out of the computational subspace, 
but we also have to take into account the geometric shift 
described in the preceding section. Indeed, since that 
shift depends on the flux drive amplitude, the resonance 
frequency of the red sideband process changes during 
the application of FC drives with G(A0(t)). In order 
to cancel that effect, we drive the sideband at the res- 
onance frequency for a geometric shift that corresponds 
to Acj)', as illustrated on Fig. |4j This amplitude is cho- 
sen such that the pulse area for which the geometric 
shift is above G(A<fi') is the same as the pulse area for 
which it is below G(A<p'). Using this definition yields 
A(f>' = Aexp(-a 2 /2cr 2 ) where 



V2a erf" 1 



-erf ( ._ 
2 \y/2a 



(72) 



Init. state 


p{\*+)) 


P(|*_» 




P(|*_» 


\4>+) 


2 x 10~ 4 


0.993 


2 x 10~ 5 


8 x 10~ 5 


10-} 


0.993 


8 x 10~ 5 


6 x 10~ 5 


8 x 10 -5 


\i>+) 


8 x 10~ 5 


9 x 10" 5 


7 x 10~ 5 


0.991 




5 x 10~ 5 


3 x 10~ 5 


0.978 


3 x 10~ 5 



TABLE II. Population transfer succes rates from the basis of 
product states to the basis of Bell states using the pulse se- 
quence displayed in Table [I] The phases are 4>i = 0.053, 02 = 
2.31, and c/> 3 = 5.59. 



by evolving the two-transmon-one-resonator system un- 
der the Hamiltonian of Eq. (42 1, along with the FC drive 
Hamiltonian for the pulse. The simulation parameters 



This method is first applied to simulate a pulse 



are indicated in Table [I] To generate the sideband pulse 
Rqi , the target qubit splitting is modulated at a fre- 
quency that lies exactly between the red sideband res- 
onance for the spectator qubit in states |0) or |1), such 
that the fidelity will be the same for both these spectator 
qubit states. We calculate the population transfer prob- 
ability for 1 1; 0) O |0; 1) after the pulse and find a success 
rate of 99.2% for both initial states |1; 0) and |0; 1). This 
is similar to the prediction from Eq. ( 41 1 , which yields 
98.7%. The agreement between the full numerics and the 
simple analytical results is remarkable, especially given 
that with \S±/e n \ = 0.23 the small S± -C e n assumption 
is not satisfied. Thus, population transfers between the 
transmon and the resonator are achievable with a good 
fidelity even in the presence of Stark shift errors coming 
from the spectator qubit (see Section III). 

The pulse sequence U cnt of Eq. (50) is then used to 
numerically prepare qubit-qubit entangled states. The 
parameters of every pulses entering in U ent are presented 
in Table [I] Given that as many as four sideband pulses 
are used, one might expect the success probability to be 
rather low. However, applying that sequence on each of 
the separable states (|00) ± |01))/V2, (|10) ± |11))/V2, 
we optain overlaps with maximally entangled states of 
~ 99%, as shown in Table [TlJ The phases coming from 
the geometric shifts, which appear in Eq. (51) describ- 
ing [Tent, are found by choosing the set {<f>i, <f>2, 4>3 } that 
optimizes the transfer rates. 

The average fidelity is found numerically by using as 
initial state the Choi matrix J2jk b')(^|e ® L?')(&|i with 
j,k G {|00;0),|01;0),|10;0),|11;0)}. The evolution is 
computed numerically on the first copy under the Hamil- 
tonian of Eq. ( 42 ) , along with the drive Hamiltonians 
for the pulses. In this way, we find the set of phases 
{4>i, 4>2, ^3 } that best fits the numerically realized gate 
and find an average fidelity of 99.1%. It should be noted 
that the presence of decoherence can only degrade these 
fidelities. In order to obtain better fidelities one would 
need to use more sophisticated pulse shaping techniques 
to compensate for the Stark shift errors. 

Finally, we included dissipation in our model using 
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FIG. 5. Average fidelity of the pulse sequence U Dn t as a func- 
tion of cavity damping damping n. The simulation includes 
realistic damping Ti = 2fj,s of the transmons and neglects 
pure dephasing. Purcell decay is taking into account in the 
simulation. The other parameters are given in Table [I] 



a Markovian master equation for two MLS and a res- 
onator Because the resonator is loaded with real (and 
not only virtual) photons, we find that the protocol is 
sensitive to photon loss k. Fig. [5] shows the average fi- 
delity as a function of the cavity decay rate. High gate 
fidelities require low cavity damping. This is possible us- 
ing 3D cavities [3 or taking advantage of the multi-mode 
structure of 2D resonators !17)]. Note that the depen- 
dence of the fidelity on qubit relaxation and dephasing is 
similar. 



VI. DISCUSSION 

Table |III| summarizes theoretical predictions and ex- 
perimental results for recent proposals for two-qubit 
gates in circuit QED. These can be divided in two broad 
classes. The first includes approaches that rely on an- 
ticrossings in the qubit-resonator or qubit-qubit spec- 
trum. They are typically very fast, since their rate is 
equal to the coupling strength involved in the anticross- 
ing. Couplings can be achieved either through direct ca- 
pacitive coupling of the qubits with strength Jq [3S], or 
through the 11-02 anticrossing in the two-transmon spec- 
trum which is mediated by the cavity [H [TU] . The latter 
technique has been successfully used with large coupling 
rates Jn_o2 and Bell-state fidelities of ~ 94%. However, 
since these gates are activated by tuning the qubits in 
and out of resonance, they have a finite on/off ratio de- 
termined by the distance between the relevant spectral 
lines. Thus, the fact that the gate is never completely 
turned off will make it very complicated to scale up to 
large numbers of qubits. Furthermore, adding qubits in 
the resonator leads to more spectral lines that also re- 
duce scalability. In that situation, turning the gates on 
and off by tuning qubit transition frequencies in and out 



of resonance without crossing these additional lines be- 
comes increasingly difficult as qubits are added in the 
resonator, an effect known as spectral crowding. 

The second class of gates includes cross-resonance, 
second-order sidebands, and our proposal. Since these 
approaches rely on external drives, their rate is not 
strictly set by the relevant coupling rate between involved 
states, but also by the ratio between the drive strength 
and some detuning, which in practice is always smaller 
than 1. This makes these gates slightly slower than those 
relying on anticrossings. Since these approaches rely on 
external drives that can be turned on and off at will, 
their on/off ratio is in principle infinite. Therefore, even 
if they take longer times than anticrossing approaches, 
they will become more interesting as gets better. Ad- 
ditionally, since they are operated either with fixed qubit 
frequencies or with qubit splittings oscillating around a 
central frequency, they should be less subject to spectral 
crowding. It should be noted that our proposal is pre- 
dicted to be faster than cross-resonance, since it scales 
with g rather than with J, which is of second order in g. 
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Appendix A: Diagonalization of the 
two-qubit-one-resonator system in the dispersive 
regime 



2 We use Eq. (2.6) of Ref. [35] with k nl 



As shown in Sec. II A after the application of the gen- 
eralized dispersive transformation of Eq. ^ , the second- 
order two-qubit plus one resonator Hamiltonian Eq. ^ is 
not completely diagonal. In this appendix, we diagonal- 
ize the remaining non-diagonal term which corresponds 
to virtual qubit-qubit interaction. 

This is done by first realizing that the Hamiltonian 
of Eq. ([5| breaks down into 4x4 blocks correspond- 
ing to a given Fock state n. One such block is given 
in Eq. |7]) and can be diagonalized exactly. Doing so, 
we get four eigenenergies and four eigenstates per pho- 
ton subspace. These eigenstates span a basis labeled 
{|00;n), |01;n), |10;n), |ll;n)} that defines two effective 
uncoupled qubits whose frequencies depend on the res- 
onator state. Expressed as a function of the eigenstates 
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Scheme 


Operation 


Theoretical rate 


Crossings 


On/off ratio 


Capacitivc coupling 36 


ViSWAP 


~ Jc 


Yes 


~ a q/ j c 


11-02 anti-crossing 181119) 


CPHASE 


~ ^11-02 


Yes 




Cross-resonance 1211 


CNOT 


~ J«/A Q 


None 


oo 


2 nd -order blue SB ITU ITSl 


Bell state gen. 


~3 3 £ 2 /A 4 


None 


oo 


1 st -order red SB 


CNOT 




None* 


oo 



TABLE III. Schemes for two-qubit operations in circuit QED. e is the strength of the drive used in the scheme, if any. (*) 
There are no crossings in that gate provided that the qubits have frequencies separated enough that they do not overlap during 
FC modulations. 



in the dispersive basis, these effective qubit states are 



1 10; n) = — sin a n \gen) + cos a n \egn), (Ala) 

|01;n) = cos a n \gen) + sm a n \egn) , (Alb) 

|00;ra) = cos (3 n \ggn) + sin/3 n |een), (Ale) 

\U;n) = - sin P n \ggn) + cos/3„|een), (Aid) 



with the mixing angles 

AS/2- Jj2 + (A«/2)2 



a,, = arctan 



j3 n = arctan 



J 



S5/2- a / J2 + (S%/2)s 
~J 



The corresponding energies are 



E 



10/01 



1= ±Jj2 + (A S /2)2 



1/00 



±^/j2 + (S™/2)2. 



(A2a) 
(A2b) 



(A3) 
(A4) 



Denning the Pauli operators t\ = (|1)(1| - |0)(0|) ® 1^ 
and t% = l' 1 * 1 ® (|1)(1| — |0)(0|), we can express the diag- 
onalized two-qubit plus one-resonator Hamiltonian as 



r « f « + ^ v /j 2 + (S Q /2 + a t a S s )2 + ^j2 + (A Q /2 + at a A s )^ ^ 
+ (^J 2 + (E Q /2 + at aSs )2 _ jj2 + (Aq/2 + a t aAs )2^ ^. 



(A5) 



In the regime where J <C £q + 2nSs, Aq + 2nAs, i.e. at 
large qubit-qubit detuning and low photon number, we 
can expand the square roots to second order in J/(Sq + 
2n£g) an d J/(Aq + 2nAg) to obtain the approximate 
Hamiltonian of Eq. ( 10 ). 



the effective qubit Pauli matrices t. 



cos 2a + cos 2/3 l cos 2a — cos 2(3 



2 sin a cos a (t}_t+ 
2 sin $ cos /3 (tAt 2 
cos 2a — cos 2/3 1 

T„ 



+/ ' 

cos 2a + cos 2/3 2 



(A6) 



Using Eq. (Ala I to (Aid), it is useful to express the 
Pauli matrices in the dispersive frame a± z in terms of 



2 sin a cos a (t_ t . 
2 sin /3 cos /3 (rir 2 



rirl), 



(A7) 
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o± = cos d cos (3 t_J_ — sin a sin (3 ri 

— cos d sin /3 r 2 — sin d cos /3 t± , 



<t_(_ = cos d cos (3 t± + sin d sin /3 tv 



cos d sin /3 rii" 2 + sin d cos /3 t_J_ -ri? , 



1 ~2 



(A8) 
(A9) 



where d and /3 are operators obtained from the mixing 
angles a n and /3 n by replacing n by a^a. 

Using these expressions, we can now express any con- 
trol drive on the qubits in the diagonal frame. As a first 
example, it is useful to consider a direct drive on qubit 1 
which is described in the bare frame by 



H dl =e dl {t) (ale™" + H.c.) 



(A10) 



Ignoring dispersive corrections, in the frame rotating at 

the qubit and resonator frequencies and assuming ujdi — 
(2) 

Wo , we find an additional term of the form 



Across = —£<ii(t) I cos d sin/3 + sin d cos /3 ) 17 r 2 . (All 



This corresponds to the cross-resonance gate [BJ [ST] . The 
present description has the advantage of including the ef- 
fects of the resonator state. Indeed, we find that the Rabi 
frequency associated to this gate depends on the photon 
number. For J <C Aq and neglecting counter-rotating 
terms, this rate is proportional to J/(AQ+2nAg). There- 
fore, depending on the sign of Ag with respect to Aq, ad- 
ditional photons act as an increased or decreased qubit- 
resonator detuning. If Ag and Aq have the same sign, 

which happens for loIP, U)a < or wi 1 \wi 2 ' 1 > w r , the 
cross resonance rate increases when the cavity is filled 
with photons. Otherwise, that rate decreases when pho- 
tons are added, allowing us to expect a saturation of the 
cross-resonance rate as a function of the drive power. 

We now express the FC drives in the diagonalized ba- 
sis. Expressing the Hamiltonian of Eq. ([I]) in the diagonal 
frame using Eq. ( A6 1 and (IA7I) , we find 



<TW = Hl{t) + Hl(t) + HUt) + H 2 SB (t) 

+ H PO (t) + H QQ (t)+H* Q (t), (A12) 

where the effects of the J coupling are taken into account 
exactly. In this expression, we have defined 



H 



^) = {4 X) / (1) (*) 
fl2(t) = {-*«/W(t) 



cos 2d + cos 2/3 



8<, 2, /< 2 >(f) 



} 



cos 2d- cos 2/3] +s! 2) / (2) 0) 



cos 2d — cos 2/3 
cos 2d + cos 2/3 I ^ 



fff B (t) = - f^(t) (A (2) cosdcos/3 + A (2) sin d sin /3) aV + (a (1) cos d cos /3 + A (1) sind sin /?) <Jn 



+ (a (2) sind cos (3 - A (2) cos d sin /§) aV+rf + (a (2) sind cos /3 - A (2) cos d sin fij aV^rf + H.c 
ffp (t) =- ( a 2 +a t2 



2 



+ 



A (i)A( 1 )/( 1 )(t) ( C os2d + cos2/3) - A (2) A (2) / (2) 0) (cos 2d - cos 2/3 

2 

A (i) A ( 1 )/( 1 )(t) (cos 2d -cos 2/3) + \W\WfW(t) (cos 2d + cos 2/3^ " 
+ 2 sind cos d (a«A« /«(<) _ A^A^ /( 2 )(t)) (rlr\ + rjr 2 ) 
+ 2 sin /3 cos /3 (a«A« /«(*) + A< 2 >A< 3 > /< 2 >(i)) (rlr 2 + r|r 2 )} , 

- /(1)( ^ /(2)( ^ 

fil)(t) \ ^^ xi cos 2/3 + (/W(M X) + / (2) W4 a) ) ^/?eos/3 



sm a cos a 



(A13) 
(A14) 



ff| B (t) = - (A (1) cosdcos^- A (1) sind sin f3\ alrl + (a (1) cos d cos /3 - A (1) sind sin $j aVj 

- (a (1) cos d sin /3 + A (1) sind cos /3j c$t\t\ - (a (1) sind cos /3 + A (1) cos d sin aV^r 2 + H.c. , (A15) 



(A16) 



(A17) 



(A18) 



"$ Q (t) =- (/ (1) (*) + / (2) (i)) {^ocosdsind (cos 2 /3t|tA + sin 2 {3tIt\} r 2 - (cos 2 /3r 2 t 2 + sin 2 /3r 2 r|) 



xi cos /3 sin /3 [(sin 2 6lt\tI + cos 2 drirj}:) r 2 - (cos 2 dr 2 r 2 + sin 2 dr 2 t 2 )] } , 



(A19) 
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where x = - A^A^ and x x = A^A^ + 

A (1) A( 2 ). We notice the presence of an additional term 



-ffqq (i) with respect to the result of Eq. ( 14 ) which repre- 
sents oscillating shifts on qubit i conditional on the state 
of qubit j. However, in the limit J <C Aq,£q, these 
terms are of fourth order in g^> and g^ and are thus 
negligible. 



Appendix B: Dispersive theory for two many-level 
systems coupled to a resonator beyond the RWA 



As explained in Section IV A the two-level dispersive 
theory of Section |H A| can be extended to the case of two 
many-level systems coupled to a resonator. The Hamilto- 



nian of this system, described in Eq. (42), can be partly 



MLS 



exp 



k) 



G 



(AO 
CR 



k=l,2 



G 



(k) 



r (k) 



fi{k) _|_ f-i{k) 
U R2 "T" Lj CR2 



(Bl) 



The first two terms in the exponential are the familiar 
terms introduced in Ref. [37], here extended to a MLS. 



= E ^ngV - nS,a (B2) 

41= E tf ) V&i*-{APy*$!i.A (B3) 



with all dispersive parameters defined in Section IV A 



Including only these terms would lead to incomplete di- 
agonalization of the Hamiltonian to second order in the 
couplings. Indeed, we would obtain a resonator squeez- 
ing term (37J and terms involving jump operators be- 
tween second-nearest-neighbour MLS states [31] . To can- 
cel these terms, we need to add second-order corrections 
to the dispersive transform. These are 



Gg?= E £fVn«-H.c., (B4) 
GS= E (efVa + ef (fe) )nW 2 -H.c., (B5) 



r (k) _ \ ~> ^'(fcWfc) 2 TT „ 



(B6) 
(B7) 



where we have defined 



(fc) 



(*) ■ (*) 



(fc) 



,(fc) 



,(A0 



■n. 



/(fc) 



2H fc) -s) 



"(fc) _ 



„(fcK(fc) 



2y 

k) 



(fc) A C=) 



2 



'7 



(fc) 



2 U 



(fc) 

+2 



t t(k) = 



CO 
'(fc) 



(fc) 



2w, 



2 w 



(fc) 

+ 2 



Co' 



(fc) 



2o>, 



diagonalized by the following dispersive transformation and 



T) 



(fc) 



Jfc)x(fc) 



'(fc) _ Jfc) A ( fe ) 



We note that 



„(fc) _ r/, ,w , ,w\ c, ,w ,, 



v (fc)\(fc) r/, ,(fc) 



(fc)> 



„(fc) x (fc) 
(fc) A (fc) 



,(fc) 



,(fc) 



(B8) 
(B9) 

(BIO) 

(Bll) 
(B12) 

(B13) 
(B14) 

(B15) 



and similarly for r^. This means that for a MLS that 
can be described as a weakly anharmonic resonator, such 
as the transmon, terms involving r\i or rj[ are negligible. 
More precisely, these terms can be safely neglected when 

\(w i+ 2 - Wj+l) - (U3 i+ l - Ui)\ < W r , V i. 

Using the Campbell-Baker-Hausdorf formula, we find 
that, to second order in the couplings, the generalized 
dispersive Hamiltonian is 



+ E 



where 
J« = \£ ] ( 



E E [^n^ + sf^Vc 

A.— 1,2 iG^ 



j -n (1) n (2) 



j/ n (l) jj(2) 



A 



A 



-,(2) 



(2) 



A 



(i) 



A 



Ha, 
(B16) 

, (B17) 
(B18) 



This Hamiltonian is not yet diagonal because of the 
Jij and Jjj terms. To complete the diagonalization, we 
assume that the detunings between transitions in the first 
and the second MLS are all much larger than and 
J'ij V i, j- Then, we can apply the unitary transformation 



exp 



A/,H (1) 



n 



(2) 



(B19) 



A J. n (i) n (2) _ H 

' ,l ij AA i,i+l AA jj + l 
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where 



~(2) 



A- 7 = 

y -(i) 



v 



i+l.i " r "'i+l.i 



(B20) 
(B21) 



To second order in the couplings, this transformation 
suppresses the interaction between the two MLS and 
brings Eq. (B16l to its diagonal form Eq. (43 1. 
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